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I. INTRODUCTION 



During the last decade a great effort has been made to understand the structure of 
matter at finite temperature and density j|. As a generally accepted hypothesis the strongly 
interacting matter at high baryon density is in a color superconducting state. Although the 
first attempts in this direction had been made already long ago Q only recently the subject 
achieved a new boost [3|. The simple observation was made that the quark-quark interaction 
at moderate densities should be of the same order as quark- ant iquark interaction. This leads 
in turn to the sizable value of the superconducting gap A ~ O(100)MeV. What is more 
interesting it was proved that because of the non-abelian structure of the underlying fields 
in QCD, the superconduction gap grows with the chemical potential at very high density 

However, the main interest is concentrated at moderate densities possible to exist in the 
core of the compact stars. Using the selfconsistent calculation in the three flavor Nambu 
- Jona-Lasinio model the phase diagram of matter was elaborated under the reasonable 
assumptions [5] . It occurs that depending on the value of the diquark coupling constant the 
2SC (or g2SC) phase occupies the large region of the phase diagram at moderate densities 
and temperatures. In the construction of the phase diagram the superconducting phases are 
usually compared to the uniform chiral phase. However, there are also other possibilities. 

It was shown that the chiral phase can develop the wave in the Lorentz scalar - pseu- 
doscalar channels which has a lower energv than the uniform state p. This effect is well 
known in the low dimensional models (eg. |7(). Nevertheless, it can be argued that the phe- 
nomenon of the chiral wave might also appear in the 3+1 dimensional QCD. The Nambu 
realization of the chiral symmetry breaking [q] is based on the idea borrowed from the BCS 
theory of superconductivity. The effective contact four-fermion interaction between quarks, 
described by the coupling constant G, is assumed to be attractive in the scalar, isoscalar 
channel (G > 0). It is then energetically favorable for the system to create the quark- 
antiquark pairs which built the Lorentz scalar chiral condensate. However QCD is invariant 
under SU{2) V x SU{2)a x U{1) symmetry. As a result the effective quark interaction also 
contains pseudoscalar, isovector channel with the same coupling constant G. If the chiral 
wave develops in both channels then, as one can show, the wave vector acts effectively as 
a constant magnetic field which couples to the spin of the quarks. Then, in the analogy 
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to the Zeeman effect, the quark energy levels split. The half of the energy branches lower 
their energy in compare to the uniform case, which depends on the quark spin projection 
on the wave vector axis. At finite density the quarks built the Fermi ladder, first filling 
the branches of the lower energy. This energy gain competes with the energy price due to 
the field gradients. There is a question of the detail calculations which contribution is more 
important. It was shown within the framework of the Nambu - Jona-Lasinio models or linear 
sigma models 0, 0] that there is a wave vector range in which the energy gain overcome 
the price. The above scenario shares similarity with the Overhauser effect [lCj which was 
already pointed out in the paper 

As one can see it is an important task to compare 2SC phases with non-uniform rather 
than uniform chiral phase. As a first att emp t the 2SC phase was compared to the non- 
uniform chiral phase at zero temperature [l2j. In this paper the calculation was extended 
to non-zero temperature. As a result the phase diagram in T — /i plane of the two flavor 
Nambu - Jona-Lasinio model including 2SC, uniform chiral and non-uniform chiral phases 
was calculated within the mean-field approximation. 



II. MEAN FIELD AND NAMBU-GORKOV FORMALISM 

We consider the NJL model at finite temperature: 

H = f {4>{il u d u + /r^V + G [{^Y + fe.rV) 2 ] + G'(^ c i l5 r 2 \ A tlj)(& l5 T 2 \ A tlj c )} (1) 

J X 

where if) is the quark field, ip c = Cip T the conjugate field and /i is the quark chemical 
potential. The color, flavor and spinor indices are suppressed. The vector f is the isospin 
vector of Pauli matrices and X A , A = 2,5,7 are three color antisymmetric SU(3) group 
generators. The integration J = J dr J d 3 x, where (3 is the inverse temperature and 
derivative operator d u = (id T , V). The coupling constant G describes the interaction in the 
isospin singlet, Lorentz scalar and the isospin triplet, Lorentz pseudoscalar, quark- ant iquark 
channels whereas G' describes the interaction in the color 3, flavor singlet, Lorentz scalar 
diquark channel. Both couplings are treated as independent. In the momentum space a 
three momentum cut-off A is introduced which finally defines the model. 

Using the Hubbard-Stratonovich transformation one can perform partial bosonization of 
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the model which leads to the standard expression 

H = - £ {i>(iYd u + yU7o + cr + ry 5 7f ■ r)ip + ?/> c (ry M <9 M - p,j + a + ■ r T )ip c (2) 
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+G'(A^75T 2 A^ + A^z 75 t 2 A a ^) - - 
We analyze the model in the mean field approximation within the following ansatz 

cr = — Mcosq-x (3) 
TT a = -M8 a3 sin q- x 
A A = A5 A2 

which takes into account three phases 

• the homogenous chiral phase Ch when the wave vector q = 0, the gap parameter A = 
and the constituent quark mass 

• The non-uniform chiral phase NCh when both the wave vector q ^ and constituent 
mass M^O whereas the gap parameter A = 0. 

• The superconducting 2SC phase when q — 0, M = and A ^0. 

The chiral phases and 2SC phase could coexist with each other, as actually is the case for 
moderate densities, which was showed at T = . One can also consider other condensates 
as for example the mixed states of 7r + , 7r~ with the zero net charge 0, 0]. All conclusions 
given below which are based on the energy considerations does not change because the model 
treats all chiral fields at the same level. 

The finite temperature analysis would be performed using Matsubara formalism in the 
Nambu-Gorkov representation. The partition function has a form 

Z-fDtD+DiJVJtolHD^D^H. (4) 

The mean field ansatz Q is space dependent thus before integration over the fermionic fields 
one makes rotation ijj' = VUip, where U = exp (275X36/* • x). Then after introduction of the 
Nambu-Gorkov (NG) basis \ T = C0'> Vc) one arrives at the mean field action 



S-\x,y) = 



(5) 



iY (d u - \vy^q v ) + //tq - M ry 5 r 2 A 2 A 

ry 5 r 2 A 2 A* vf{d v - \^l^qv) - jijo - M 



where q u = (0, q). The bilinear form of Hmf can be integrated over the fermionic fields with 
the mean field result 



Z M f = exp 



, / M 2 |A| 2 \ 1 l 



Determinant of the inverse propagator is calculated in the momentum space 

S-\K) = [ S-\x,y)exp(-iK(x-y)), 
Jk 



where 



T E 



d 3 k 



, K v = (—iuj n , k), uo n = (2n + 1)ttT. 



(6) 



(7) 



(8) 



* k — J (2tt)- 

n=— oo v ' 

The inverse propagator S~ X {K) for any given K is a 48 x 48 matrix in the color, flavor, Dirac 
and NG space. However, it decays into two block-diagonal sub-matrices: 16 x 16 matrix 
S^and 32 x 32 matrix S^ 1 . This simplification follows because the diquark interaction 
does not mixed the blue and red/green quarks. The calculation of these two determinants 
is straightforward although rather tedious and leads to the result 



det S^ 1 



n 1 

.i,k=± 



f 2 - K 2 ' 



(9) 



e+,± = V(// + £±) 2 + |A| 2 , e_,± = £±) 



IAI 



<f 2 



A; 2 + M 2 + ^- ± y (<f • k) 2 + M 2 cp 

where det Sq 1 = det[5^_ ]. After substitution of (JHJ) into (0) one gets the grand potential 
Q = —T In Zm f in the form 
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v 

Vo 



V -AT 
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'1 + exp — ^ ) H — In -(1 + exp — - -J 
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(10) 



M 2 IAI 2 ^ f A d 3 k , . ^ f d 3 k /rn 



4G 4G 
e° fc = e lk (A = 0). 



In the limit of zero temperature the second term of the potential Q vanishes. It means that 
the term Vq describes the zero temperature contribution. This term can be rewritten in the 
form useful for further calculation 



To 



M 2 |A| 2 M 2 F 2 q 2 
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FIG. 1: The phase diagram in the T — \x plain in GeV for the coupling constants G' = G/2 (left 
panel) and G' = 3G/4 (right panel). 



where the expansion in the power of the wave vector q was used as described in |12[ and 
M = 0.301 GeV is a constituent quark mass at zero density. The quantity E = Ei(q = 0) 
and F,r = 93 MeV is a pion decay constant. 

III. RESULTS 

The phase diagram describes the global minima of the potential Q as a function of tem- 
perature T and chemical potential /i. Further constraints can be also imposed on the system 
as color and charge neutrality. We do not include them in this paper because its influence on 
the competition between chiral and 2SC phases is not quantitatively substantial. This prob- 
lem would be discussed at the end of this section. The thermodynamic potential Q = Vq + Vt 
is minimized with respect to M, q and A. The zero temperature contribution Vq is given by 
equation (fTTjl and Vr is the finite temperature part of expression (fTTH) (second term of fl). 



The model parameters can be fitted to the values of the chira 



condensate and the pion mass 



which give C ,-■ r,.()l CVV-- and the cut-off A = 0.65 GoV Q- 

The coupling constant G' 

was treated as independent from G. The minima were found using standard Nelder-Mead 
Simplex algorithm. The calculation was performed starting from three different points re- 
lated to three different phases. The global minimum was that of the lowest value of the 
thermodynamic potential. 

In the Fig. 1 the phase diagram for G' = G/2 (left panel) and G' = 3G/A (right panel) 
are described. There are five different phases 
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• the quark matter phase (QM) with M = 0, q = 0, A = 0. 

• The chiral phase (Ch) M ^ g = 0, A = 0. 

• The non-uniform chiral phase (NCh) M ^ 0, q ^ 0, A = 0. 

• The mixed phase of color superconducting and non-uniform chiral phases (SNCh) 
M^O, gVO, A^0. 

• The color supercoducting phase (2SC) M = 0, q = 0, A ^ 0. 

Each of these phases differ by the order parameter and can be easily distinguished 
from each other. The phase transition lines between Ch/NCh, QM/NCh, Ch/SNCh and 
2SC/SNCh are first order. This is connected to the fact that the entrance to the region of 
non-uniform condensate is associated with the jump in the value of the wave vector. This 
discontinues change entails the discontinues changes in the mass and gap parameters (see 
Fig. 2). The QM/Ch, QM/2SC and NCh/SNCh lines are second order phase transitions. 
This is nicely seen on the right panel of the Fig. 2 which depicts the values of order parame- 
ters at the temperature T = 0.035 GeV as a function of the chemical potential. First there 
is a 1st order phase transition from Ch into the NCh phase around fx = 0.29 GeV which 
is followed by the continuous phase transition into SNCh phase around \x = 0.33 GeV (rise 
of the gap parameter A). Eventually there is the 1st order transition into the 2SC phase 
at /i — 0.35 GeV, where both mass and the wave vector parameters vanish. The change of 
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FIG. 2: The value of the mass M (stars), wave vector q (boxes) and the gap parameter A (triangles) 
for T = 0.075 GeV (left panel) and T = 0.035 GeV (right panel) for coupling constants G' = 3G/4. 
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the character of the QM/Ch transition line from the second to the first order takes place in 
the merge point around (/i, T) = (0.221, 0.125) GeV where QM, Ch and NCh phases meets 
with each other. 

The details of the phase diagram depends on the strength of the diquark coupling con- 
stant. The lower value of G' the smaller part of the diagram is covered with the pure 2SC 
phase. In the Fig. 1 it is seen that at low temperature the SNCh/2SC phase boundary starts 
at fi m 0.409 MeV for G' = G/2 in compare to y. » 0.35 MeV for G' = 3G/4. However when 
temperature grows both phase diagrams start to resemble each other (also quantitatively). 
Such a behavior is expected because the color superconductivity is rather low temperature 
phenomenon and does not influence the phase structure at higher temperature. In particular 
the region of NCh field is almost independent of the value of G' coupling constant. The only 
difference is that SNCh phase extends to the wider region. Also the melting temperature 
of the 2SC phase is larger at the right panel compare to the left one. This follows from the 
fact that the gap parameter takes larger values when G' = 3G/4. 

The phase diagram at higher temperature (T « 120 MeV) has to be treated with caution 
because of the corrections from the pion loops. In particular the characteristic insert of the 
NCh phase (the horn going to the left from the merge point of NCh/QM lines) possesses the 
wave vector of a small value \q\ -C irT thus the temperature fluctuations rather destroyed 
that part of the phase. Moreover the details of the phase diagram at higher temperature 
certainly depends on the regularization procedure, as one can conclude comparing Fig. 
1 with similar from the paper 9}. The difference between these two approaches comes 
from the different way of the regularization of the Vr part of the thermodynamic potential. 
Nevertheless at the lower temperature the results are robust within the expected accuracy 
of the NJL model. 

Finally one has to remember that 2SC phase describes the colorful strongly interacting 
matter. The color charge neutrality condition can be easily accommodated through the 
additional color chemical potential (iJ). The 2SC phase is suppressed under this condition 
thus as a result the (S)NCh phases extend to the larger region of the diagram. However the 
change in the phase diagram structure is not substantial. At zero temperature it differs only 
in details at the few per cent level at most Q|. 

The influence of the beta equilibrium condition on the superconducting phases can result 
in the change of the phase structure of the color superconductor ^| . Then one needs rather 
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to compare the non-uniform chiral phase with the LOFF superconducting phase. However, 
this task certainly requires different study and is outside the scope of the present paper. 



IV. CONCLUSIONS AND FUTURE PROSPECTS 

In this paper the non-uniform chiral condensate NCh was compared to the uniform 2SC 
phase (j3J) at finite densities and temperatures. It was shown that the NCh phase exists in 
the temperature range between 20 - 100 MeV above chemical potential /i ~ 250 MeV. It 
also coexists with 2SC phase at lower temperatures. The range of the phase coexistence 
depends on the relative strength of the chiral and diquark couplings. For G' = G/2 the 
density window open at /z ~ 0.3 GeV and lasts for about 100 MeV whereas for G' = 3G/4 
it shrinks to about 50 MeV. The detailed phase diagram can be found in the Fig. 1. It 
is worth to remind that NCh phase can have an important influence on the neutron star 
structure and also can be a source of the strong magnetic field js ( 12|- Then, the mixed state 
SNCh can have very unusual features as composed of NCh and 2SC phases. This interesting 
expectation remains to be revealed in some future work. 

The light quarks are assumed to be massless in the model. The introduction of the non- 
zero mass leads to the tilt in the chiral wave as was shown in the low dimensional model 
0. There is not known exact solution to this problem in 3+1 dimension. However, the 
effect has to be proportional to the current quark masses, which are the smallest energy 
scales in the system and should not have an important influence on energetic considerations. 
Nevertheless, the detailed analysis remains as an open question. 

Another important problem is the comparison of the NCh phase and the color supercon- 
ducting LOFF phase. This last phase appears naturally under the condition that matter is 
at beta equilibrium in the interior of the compact stars jl^ . It is also consider as a remedy 
for the chromomagnetic instabilities in two flavor color superconductor . At the first step 
one can try the simplified version of the problem when only two wave vectors for chiral and 
superconducting phase are assumed. The more general ansatz would be more ambitious, 
however, much harder to solve. Nevertheless the full analysis of the ground state at high 
baryon densities would require the study of the non-uniform phases of the chiral as well as 
diquark condensates. 
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